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We present a phenomenological model for 2D turbulence in which the energy spectrum obeys 
a nonlinear fourth-order differential equation. This equation respects the scaling properties of the 
original Navier-Stokes equations and it has both the —5/3 inverse-cascade and the —3 direct-cascade 
spectra. In addition,our model has Raleigh- Jeans thermodynamic distributions, as exact steady state 
solutions. We use the model to derive a relation between the direct-cascade and the inverse-cascade 
Kolmogorov constants which is in a good qualitative agreement with the laboratory and numerical 
experiments. We discuss a steady state solution where both the enstrophy and the energy cascades 
are present simultaneously and we discuss it in context of the Nastrom-Gage spectrum observed 
in atmospheric turbulence. We also consider the effect of the bottom friction onto the cascade 
solutions, and show that it leads to an additional decrease and finite-wavenumber cutoffs of the 
respective cascade spectra which agrees with existing experimental and numerical results. 

PACS numbers: 



I. Leith differential model and its relatives 

In 1967 paper 0, Leith proposed the following differ- 
ential equation model as a model of 3D turbulence, 

where t is time, k is the absolute value of the wavenum- 
ber, v is the kinematic viscosity and the one-dimensional 
energy spectrum E(k, t) is normalised so that the kinetic 
energy density is J E dk. The third term on the RHS 
of Eq. Q is obvious and describes the linear process of 
viscous dissipation. The first term contains a nontriv- 
ial model of the nonlinear processes in turbulence which 
rests on three basic properties: 

1. The total energy density / E dk is conserved by the 
inviscid dynamics. The characteristic time of the spectral 
energy redistribution for local interaction of turbulent 
scales is of order of the vortex turnover time, l/Vk 3 E. 

2. The steady state in forced turbulence corresponds 
to a constant energy cascade though the inertial range of 
scales which is described by the Kolmogorov spectrum, 

E = CP 2 / 3 k- 5 / 3 , (2) 

where P is the energy flux (constant in t and k) and 
C « 1.6 is the Kolmogorov constant. 

3. When the wave- number range is truncated to a 
finite value and both forcing and dissipation are absent, 
turbulence reaches a thermodynamic equilibrium state 
characterized by equipartition of energy over the wave- 
number space |2| ■ In terms of the one-dimensional energy 
spectrum this means E oc k 2 which is obviously a steady 
state solution of the equation (JTJ for v = 0. 

Even simpler first order differential equation was pro- 
posed by Kovazhny, which can model the cascade prop- 
erties but ignores the thermodynamic distributions. Due 



to their simplicity, Leith and Kovazhny models and their 
generalisations to include other physical effects turned 
out to be very efficient. In particular, they were used for 
description mixed cascade-thermo solutions correspond- 
ing to the bottleneck effect 0], superfluid turbulence in 
3 He with friction against the normal component [5(, in- 
fluence of heavy particles on turbulent suspensions jj| . 

Differential models have also been proposed and used 
before for wave turbulence. In particular, Hasselmann 
and Hasselmann Q and, independently, Iroshnkov Q 
(thereafter HHI) proposed a model for the spectrum of 
deep-water surface gravity waves, which is a fourth-order 
nonlinear differential equation. The order of this equa- 
tion is higher than in the Leith model because, in addi- 
tion to the direct energy cascade and the thermodynamic 
energy equipartition, it describes the inverse cascade and 
the equipartition state of waveaction, - an additional in- 
variant of the water wave system. Recently this model 
was used to describe a sandpile-like behaviour energy cas- 
cades of waves in discrete A:-space . Similar model was 
proposed for turbulence of Kelvin waves on thin quantum 
vortices, and it was used to study sound radiation by such 
waves |9( • This model is also of the fourth order because 
the Kelvin wave system also conserves waveaction |30j . 



II. Differential Model for 2D turbulence 

A 2nd-order differential model for 2D turbulence was 
proposed by Lilly 0], see equation 1|11JI below. It has 
both the energy and the enstrophy cascade solutions, but 
it ignores the thermodynamic equilibria. As we will see 
later, this causes a significant defect in description be- 
cause it fails to predict the important difference in val- 
ues of the Kolmogorov constants for the direct and the 
inverse cascade spectra. The minimal model that would 
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include both cascades and the thermodynamic solutions 
has to be of the fourth order. In some sense, such model 
should be a hybrid of the Leith model, since it has to re- 
spect the Navier-Stokes scalings, and the HHI model, be- 
cause of presence of the additional invariant - enstrophy. 
We therefore postulate the following 4th-order evolution 
equation for the energy spectrum in 2D turbulence, 
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where x — k 2 and D is the dimension of the system. 
Although our main object is 2D turbulence, D — 2, we 
would like to keep D as a parameter for the purposed 
which will be clear later. For v = this equation con- 
serves the mean energy density, 



(v 2 ) 



Edk 



Edx 



(4a) 



and the mean enstrophy density, 



i((V x v) 2 ) = J k 2 Edk = ~f y^Edx . (4b) 

Equation can be written as a continuity equations for 
the energy, E, and the enstrophy, II = k 2 E, spectra: 

dE dP dP dQ 

where P and Q are the flux of energy and the flux of 
enstrophy respectively, 
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III. Steady state solutions 

Equation has a general thermodynamic Raleigh- 
Jeans solution in which both the energy flux and the 
enstropy flux are zero. It is determined by the condition 
R = which yields, 



E=(Tk u - 1 )/(k 2 + ti), 



(7) 



where constants T and /x are temperature and chemical 
potential respectively. 

In addition, Eq. © has two cascade solutions: an in- 
verse energy cascade solution (P = const.), 

E = C P (-P) 2 / 3 AT 5 / 3 , (8a) 

and a forward enstrophy cascade state (Q — const.), 

E = C Q Q 2 ^k- 3 , (8b) 



where Cp and Cq are Kolmogorov constants. 

Substituting Eqs. JBJl into Eqs. JJjJl, one can find rela- 
tions between Cp, Cq and a: 

72 = a C P (3D + 2) (3D - 4) , (9a) 
8 = a C Q D(D + 2) . (9b) 

This implies the following relation between the Kol- 
mogorov constants, 



Cq 
C P 



(3D + 2)(3D-4) 



-i 2/3 



9 D(D + 2) 



(9c) 



Note that this expression turns into zero for D = 4/3 
which means that the the for this dimension the in- 
verse cascade rate turns into zero and the correspond- 
ing inverse-cascade spectrum degenerates into thermo- 
dynamic energy equipartition. Correspondingly, one can 
expect that for D = 4/3 Gaussian statistics for this range 
of scales (up-scale from the forcing scale) . 
For D — 2 we have 



Cq 
C P 



2/3 



0.367. 



(10) 



Note that D = 2 is rather close to the critical dimen- 
sion D = 4/3 and, therefore, statistics of turbulence in 
the inverse-cascade range should be expected to be quite 
close to Gaussian. This observation was previously made 
in 10] . Another consequence of the fact that D = 2 is 
close to the critical dimension D = 4/3 is that Cq is 
significantly less than Cp as seen in l|10[l. Experimental 
values for C Q are 1.4±0.3 [13 where as DNS give 1.6±0.1 
[13]] and more rece ntly 1.9±0.1 |l4|. For Cp, experiments 
give Cp « 6.5 ± 1 |l2( which is consistent with DNS (e.g. 
Cp w 7 in 0)- More recently it was argued however 
[l5j that the energy condensation at largest scales makes 
it difficult to measure Cp and its fit varies from 4.5 near 
the forcing scale to about 7.2 near the peak of the spec- 
trum. Taking a ratio corresponding to most recent high- 
resolution DNS we have C Q /C P » 1.9/6 ps 0.32 which 
is in a godd qualitative agreement with (fl"U|) . Note that 
although there is still a significant spread in the available 
experimenta and numerical data on the Kolmogorov con- 
stants, their ratio is always significantly less than one and 
our model predicts and explains this behavior. 



IV. Two-flux solutions 

Pure single cascade states © are mathematical ideal- 
isations corresponding to infinitely wide inertial ranges, 
whereas in both numerical simulations and nature these 
ranges are finite and typically both the energy and the 
enstrophy cascades are present in important turbulent 
scales. In particular, such two-cascade states are often 
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applied for explaining the Nastrom-Gage spectrum ob- 
served in atmospheric turbulence: with —3 exponent at 
low wavenumbers and —5/3 at the high wavenumber end 
(T7L Hil ITg| . To find a analytical expression for such 
two-flux spectra, Lilly introduced a second-order equa- 
tion model 1!) which ignores the thermodynamic states 
but describes both the inverse and the forward cascade 
states. 



2E - h— 

dt ~ dk 
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where b is a dimensionless constant which has to be found 
from an experimental or numerical value of Cp (using Cq 
is less desirable due to certain lack of precision in finding 
this constant related to nonlocality and respective log- 
correction of the direct cascade state). We have, 
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0.05 . 



(12) 



Equation (|ll|l is written in a form on the energy conti- 
nuity equation. It can also be re-written as an enstrophy 
continuity equation, 
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Such a 2nd-order model is less realistic than the 4th-order 
model because it gives Cp — Cq due to absence of ther- 
modynamic states. As we will see below, this also makes 
it less realistic for description of the two-flux states. 
The general steady-state solution of 1)11(1 is 



E = fc" 
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(14) 



This spectrum exhibits behaviour characteristic of the 
Nastrom-Gage spectrum: —3 exponent at the left and 
—5/3 right sides of the inertial range 0,0]. Formally, 
Cp = Cq, and substituting this equality into l(Tljl we 
would recover a solution originally obtained by Lilly [l9| 
(whose main goal was to explain the Nastrom-Gage spec- 
trum). However, we emphasize that the more precise 
4th-order model, as well as the DNS and experimental 
data, give significantly different values for Cp and Cq, 
and therefore we can expect that Eq. (fHjl to work better 
with Cp and Cq taken from the 4th-order model, or from 
the available numerical or experimental data discussed 
above. Transition between the two exponents occurs at 
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As it was pointed out in [l9j , this transition occurs with- 
out any sinks of turbulence in between of the two sources. 

Formula l|15|l gives correct asymptotics for P ^> Q and 
for P <C Q- However, by substituting it into more pre- 
cise fourth-order model we see that both the energy and 
the enstropy fluxes ((6a(l experience order of magnitude 



variation in the region P ~ Q instead of remaining ap- 
proximately constant as they should. Thus, we conclude 
that (|15fl does not work well in the transitional region 
P ~ Q and a better description should be expected from 
finding the two-flux states directly in the 4th-order model 
©. Using Eqs. ©, we can write the following differential 
equation for the general two-flux steady state, 



VES X ^/ 2 |1 (^) = -2{Px + Q)/a 



(16) 



This equation cannot be solved analytically, but it can 
be easily integrated numerically for any particular values 
of P and Q and compared to observational or DNS data. 

Note, however, that the two-source scenario is not the 
only suggested explaination of Nastrom-Gage spectrum 
and, in particular, the steeper slope at small wavenum- 
bers can appear due to non-universality because of inefi- 
ciency of large-scale dissipation resulting in condensation 
of energy at the large scales (see detailed discussion and 
references in review Hi]]). 



V. Effect of friction 

In many situations, particularly atmospheric turbu- 
lence and laboratory experiments, the bottom friction 
is believed to have an important effect on 2D turbulence, 
and in this section we will address this issue using our dif- 
ferential model. Our interest here is partially motivated 
by a question if the bottom friction can help to form a 
spectrum of the Nastrom-Gage shape, particularly in the 
view of results of || where a spectrum with transition 
from the —3 to the —5/3 exponent was obtained in 3D 
superfluid turbulence with friction (and see comments at 
|31|). For all that follows one could use the fourth order 
equation, but we chose to work with the second order 
model here because of easier algebra. Adding friction 
and ignoring viscosity in ((11(1 we have, 
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where 7 is a fc-independent friction frequency. One could 
try a power law substitution, E = Ck v , for which the 
power counting immediately gives y = —3 |23| . but the 
constant C turns out to be undefined (infinite). This is 
natural because the —3 shape coincides with the direct 
cascade spectrum in the absence of friction and, thus, the 
nonlinear term in 1(17(1 turns into zero. However, there 
should be a corrected solution which corresponds to a di- 
rect enstrophy cascade in presence of friction. Close to 
the forcing scale fc^ this correction can be found substi- 
tuting E = C Q Q 2/3 fc~ 3 (l + S Q ) into ((HI) which gives 
after linearisation, 



Z r 1/2 < 
2b C Q 



f/ 3 ln(fc/fc F ) 



(18) 
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One can see that this correction is negative and growing 
in magnitude along the cascade to high k's. Far from the 
forcing scale, fc 3> kp we gain a logarithmic factor, 
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(19) 



where /cq is a constant wavenumber at which an abrupt 
cutoff of the spectrum occurs. Note that the finite cut- 
off is a result of using the differential model and it indi- 
cates that in more general models the spectrum, although 
always non-zero, decays at these scales faster than any 
power of k. The value of &q can be obtained via numer- 
ical solution for the full spectrum matching the regions 
close and far from the forcing scale. Note that because 
neither 7 nor Q contain the length dimension the expres- 
sion for kQ must contain the forcing wavenumber - the 
only relevant quantity containing the length dimension. 

Result l|19|l about the log-correction is natural keep- 
ing in mind the degeneracy of coincidence of the direct 
cascade exponent, —3, with the exponent arising from di- 
mensional balancing of the friction and the inertial terms 
in the energy balance equation. However, papers [islE^ 
present a numerical and experimental evidence that the 
friction effect on the direct cascade is to modify the spec- 
trum exponent rather that to log-correct it. Plausible 
explanations for the descrepancy with the prediction of 
the differential model could be: (i) energy condensation 
and nonlocal interactions with large scales may have been 
important in [3, (they are not described by the dif- 
ferential model which is super-local) or (ii) insufficient 
fitting interval (one decade) in |28l . did not allow to 
distinguish between the log and power corrections. Thus, 
additional numerical experiments at higher resolution are 
desirable to obtain a larger fitting interval. 

Summarising results of (|18f) and (|19f) . we conclude that 
friction leads to a faster decay of the inverse cascade spec- 
trum and its eventual abrupt cutoff. This means in par- 
ticular that it would never take a flatter -5/3 slope at 
high k characteristic for Nastrom-Gage spectrum. 

Let us now find the friction correction to the inverse 
cascade spectrum close to the forcing scale by substitut- 
ing E = C P P 2 / 3 k- 5 / 3 (l + S P ) into Q [with added 
friction term as in (|17Jl ]. After linearisation, we have 



^ = _^ C V2pl/3 A .-2/3 



(20) 



One can see that this correction is negative and growing 
in magnitude along the cascade to low fc's. Easy to see 
by inspection that the inverse cascade spectrum also has 
a finite cutoff close to which we have 



kp) 2 . 



(21) 



Thus we see that in the inverse cascade range one could 
never have slope —3 at low k characteristic for Nastrom- 
Gage spectrum. Exact position of the cutoff can be found 



only by matching the large-scale and the small-scale parts 
of the solution, which can be done numerically. However, 
up to an order-one factor, it can be uniquely found from 
the dimensional argument which gives, 



A 3/a/|p|i/ 2i 



(22) 



In fact, existence of such an abrupt cutoff was suggested 
by Lilly [2^] that the Kolmogorov relation between the 
spectrum and the energy flux persist. This assumption 
is valid only when friction affects the flux only weakly, 
and therefore it is natural that Lilly's expression for the 
spectrum agrees with ours near the forcing scale (|2U|) but 
it fails to predict correct behaviour near the cutoff scale 
(|21|l . Estimate l|22|l was also obtained by Manin who con- 
sidered stability of large-scale vortices j^. This cutoff 
scale was also discussed in experimental works [24I Esl Efij ] 
as well as in numerical works ([2lJ and references therein) 
in the context of arrest of the inverse energy cascade by 
bottom friction. 

Combining the results for the direct and the inverse 
cascade spectra, we conclude that when forcing is present 
at only one scale the Nastrom-Gage spectrum could not 
be explained by the friction effect neither in the direct 
cascade nor in the inverse cascade ranges. Thus, 2D 
turbulence with friction is fundamentally different from 
the 3D case, and the origin of this difference can be 
attributed to presence of the enstrophy cascade (with 
exponent coinciding with the one obtained by the formal 
power counting in the in the frictional system) and to 
the fact that the energy cascades in 2D and 3D are in 
the opposite directions with respect to each other. 

Summary. In this Letter, we presented a differen- 
tial models for 2D turbulence based on the 4th-order 
Eq. ©. Based on Eq. ® we derived the ratio of the 
Kolmogorov constants 0.37 which qualitatively agrees 
with an experimentally and numerically measured value 
of about 0.32. Of particular importance is the fact that 
the model predicts the inverse-cascade constant to be 
significantly greater than the direct-cascade constant. 
Based on a reduced 2nd-order model l|ll|) . we discussed 
a mixed two-cascade solution ()14f) which was previously 
suggested by Lilly for explaining the the observed 
Nastrom-Gage spectrum of atmospheric turbulence. 
We showed that, inspite of giving correct asymptotic 
expressions in the small and large k limits, this solution 
is inacurate in the intermediate range. We concluded 
that the 4th-order model should be used instead for 
finding the two-cascade states, and doing so leads to 
the 2nd-order nonlinear ODE (|16l) which can be solved 
numerically in each particular case with given values of 
P and Q. We also examined the effect of friction on 
the inverse and direct cascades, and showed that in the 
case with only one forcing region present, the friction 
cannot lead to the Nastrom-Gage shape. Instead, for 
both cascades friction acts to reduce the spectra and to 
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terminate them abruptly at finite cutoff wavenumbers, 
which agrees with previous numerical and experimental 
studies. However, our result that friction acts to intro- 
duce a log-correction on the direct cascade rather that 
to modify its exponent is different from the conclusions 
of previous studies [2^, |2j|. A possible explanation of 
this descrepance is that non-local effects were important 
in simulations and experiments of [isl . 

Acknowledgements. This work was done as the re- 
sult of our collaboration during the the Warwick Sympo- 
sium workshop "Environmental Turbulence: from Clouds 
through the Ocean" at Coventry, UK on March 2006, and 
the ULTI Users Meeting at Lammi, Finland on April 
2006. We are acknowledge the support of this project 
by Large Scale Installation Program ULTI of the Euro- 
pean Union (contract number: RITA-CT-2003-505313), 
EPSRC Turbulence Symposium grant, the European Sci- 
ence Foundation Program COSLAB and the US-Israel 
Binational Science Foundation. 



[1] C. Leith, Phys. Fluids 10, 1409 (1967); Phys. Fluids 11, 
1612 (1968). 

[2] Two-dimensional turbulence, R. Kraichnan and D. Mont- 
gomery, Reports on Progress in Physics, 43, 547 (1980). 

[3] C. Connaughton and S. Nazarenko, Warm Cascades and 
Anomlous Scaling in a Diffusion Model of Turbulence, 
Phys. Rev. Lett. 92, 044501 (2004). 

[4] V.S. L'vov, G. Ooms and A. Pomyalov, Effect of particle 
inertia on the turbulence in suspensions, Phys. Rev. E, 
67, 046314 (2003). 

[5] V.S.Lvov, S.V. Nazarenko and G.Volovik, Energy spectra 
of developed superfluid turbulence, JETP Letters 80 (7) 
pp 535-539 (2004). 

[6] S. Hasselmann and K. Hasselmann, Computations and 
parametrizations of the nonlinear energy transfer in grav- 
ity wave spectrum. Part 1, J. Phys. Oceanogr., 15, 1369- 
1377 (1985). 

[7] R.S. Iroshnikov, Possibility of a non-isotropic spectrum 
of wind waves by their weak interaction, Soviet Phys. 
Dokl, 30, 126-128 (1985). 
[8] S.V. Nazarenko, Sandpile behavior in water wave tur- 
bulence, J. Stat. Mech. (2006) L0200 2 doi:10.108 8/1742- 
5468/2006/02/L02002 (arXiv: nlin.CD/0510054 ). 
[9] S.V. Nazarenko, Differential approximation for Kelvin- 
wave turbulence, JETP Letters, 83, pp 198-200 (2005) 
(arXiv: cond-mat/0511136 ). 

[10] V.S. L'vov, A. Pomyalov and I. Procaccia, Quasi- 
Gaussian statistics of hydrodynamic turbulence in (| +e) 
dimensions, Phys. Rev. Lett., 89, 064501 (2002). 

[11] J. Paret, M.-C. Jullien and P. Tabeling, Vorticity statis- 
tics in the two-dimensional enstrophy cascade, arxiv: 
physics/9904044 vl 21 Apr 1999. 

[12] J. Paret and P. Tabeling, Experimental observation of 
the two-dimensional inverse energy cascade, Phys. Rev. 
Lett., 79 (21) 4162 (1997). 



[13] V. Borue, Phys. Rev. Lett., 71, 24, 3967 (1993). 

[14] T. Ishihara and Y. Kaneda, Energy spectrum in the en- 
strophy transfer range of two-dimensional forced turbu- 
lence, Phys. Fluids 13 (2) 544 (2001). 

[15] S.Danilov and D.Gurarie, Nonuniversal features of forced 
two-dimensional turbulence in the energy range, Phys. 
Rev. E 63, 020203(R) (2001). 

[16] L.M. Smith and V. Yakhot, Bose condensation and small- 
scale structure generation in a random force driven 2D 
turbulence, Phys. Rev. Lett., 71 352 (1993). 

[17] G.D. Nastrom and K.S. Gage, A climatology of atmo- 
spheric wavenumber spectra of wind and temperature 
observed by commercial aircraft, J. Atmos. Sci. 42 pp 
950-960, (1984). 

[18] G.D. Nastrom, K.S. Gage and W.H. Jasperson, The at- 
mospheric kinetic energy spectrum, 10° — 10 4 km, Nature 
310 pp 36-38 (1984). 

[19] D. K. Lilly, Two-Dimensional Turbulence Generated by 
Energy Sources at Two Scales Journal of the Atmospheric 
Sciences: Vol. 46, No. 13, pp. 2026-2030 (1989). 

[20] Maltrud, M. and G.K. Vallis, Energy spectra and coher- 
ent structures in forced two-dimensional and geostrophic 
turbulence. J. Fluid Mech., 228, 321-342 (1991). 

[21] S.Danilov and D.Gurarie, Quasi-two-dimensional turbu- 
lence, Physics- Uspekhi 43 (9) 863-900 (2000). 

[22] D. K. Lilly, Numerical simulation studies of two- 
dimensional turbulence: Part 1. Models of statistically 
steady turbulence., Geophys. Fluid. Dyn. 3 289 (1972). 

[23] L.M. Smith and V. Yakhot, Finite Size Effects in Forced 
Two-Dimensional Turbulence, J. Fluid. Mech., 274 115 
(1994). 

[24] Colin de Verdiere, A. Geophys. Astrophys. Fluid Dyn. 
15, 213-251 (1980). 

[25] J. Sommeria, Experimental study of the two-dimensional 
inverse energy cascade in a square box J. Fluid Mech., 
170, p.139-168 (1986). 

[26] Dolzhanskii, F.V., Krymov, V.A. and Manin, D.Yu., 
1990. Stability and coherent structures in quasi-two- 
dimensional shear flows. Sov. Phys. Uspekhi, 33 (7), 495- 
520. 

[27] Manin, D.Yu., 1990. Characteristic Size of Vortices in 
Developed Quasi-Two-Dimensional Flows. Izvestiia. At- 
mospheric and oceanic physics., 26 (6), 426. 

[28] G. Boffetta, A. Cenedese, S. Espa and S. Musacchio, Ef- 
fects of friction on 2D turbulence: An experimental study 
of the direct cascade, Europhys. Lett., 71 (4), pp. 590- 
596 (2005). 

[29] L. Biferale, M. Cencini, A. Lanotte and D. Vergni, In- 
verse velocity statistics in two dimensional turbulence, 
Phys. Fluids 15 1012, (2003). 

[30] Waveaction is an adiabatic invariant (approximately) 
conserved by weakly nonlinear systems whose resonant 
interactions are of even order. For example, on the deep 
water the leading resonances are four-wave whereas for 
the Kelvin waves they are six-wave. 

[31] Similar question was examined by Elef Gkioulekas in his 
talk at the Warwick Symposium workshop " Universal 
features in turbulence: from quantum to cosmological 
scales" (December 5-10, 2005). However he used a 3D 
differential model which is less relevant for this case than 
the 2D model introduced in the present paper and there- 
fore we re-examine this question using our model. 



